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We introduce a new paradigm of (electromagnetic) 2D ray-chaos, featuring both guided and scat- 
tered rays in a dielectric layer with exponentially tapered refraction index backed by an undulated 
conductive surface, and illustrate its relevant features. Numerical simulations of the corresponding 
full-wave solution indicate that the system complies with Berry's conjecture in the asymptotic short 
■ \ - wavelength limit. 

(N ; 

PACS numbers: 05.45.Mt, 42.15.Dp, 42. 25. Fx 

CJO! I. INTRODUCTION 

3 ■ 

Ray-chaos (exponential separation of nearby-originating rays, due to the geometrical features of the reflecting 
boundaries, and/or the refractive properties of the medium) is a lively subject, as witnessed by an increasing number 
C\| ■ of topical papers. 

Well known paradigms of two-dimensional (henceforth 2D) ray-chaos include the Sinai |l| and Bunimovic Q billiards 

Q (internal boundary value problems) and the n-disk (n > 3) pinball-scatterer (external boundary value problem) [0. 
These systems have been extensively studied both theoretically [5J, in view of their relevance in the context of 
>mJ , quantum chaology B, and experimentally (see e.g. ||, and 0), thanks to the formal equivalence between the 2D 
Schrodinger equation for a point particle in an infinite-wall potential, and the 2D (scalar) Helmholtz equation with 
perfectly conducting boundaries of Electromagnetics, which makes an experimental approach based on microwave 
1 measurements viable. 

In this paper we introduce and analyze a new intriguing paradigm in ray-chaotic electromagnetic (henceforth EM) 
■ boundary-value problems. In its simplest (2D) form, the problem's geometry is sketched in fig. 1, and consists of a 
dielectric layer with exponentially tapered refraction index matched at z — to the outer vacuum (viz., n(0) = 1), 
[• — , and backed by a smooth perfectly-conducting periodic (undulated) surface z = Q(x). 

This geometry bears several elements of novelty. It is, perhaps, the simplest (2D) internal/external electromagnetic 
ray-chaotic boundary value problem featuring both guided and scattered rays at the same time. It also includes a 
(linear, isotropic, time- invariant, non-dispersive) dielectric layer as a basic ingredient at variance of billiards and 
pinballs, which are pure perfect-conductor-and-vacuum problems. 

The 2D ray trajectories bear some similarity to the space-time ones of the well known (one-dimensional) Fermi- 
— ' Pustil'nikov billiard In this analogy, the gradient of the refractive-index plays the role of gravity, bending the 
ray-paths downward in the dielectric layer, while the undulated surface plays the role of the vibrating table, changing 
the slope of the ray at each bounce. 
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II. RAY TRACING 

The system in Fig. 1 can be parametrized in terms of the dimensionless variables nf = n(—h), a/h, A/a. In this 
section we illustrate the main relevant results of straightforward ray-tracing simulations run for the geometry of Fig. 
1, with a/h — 1 and A/a — 2n/5. 

In Fig. 2 three examples are shown of (multi-hop) ray-paths originating from rays with close-by incidence points 
(Ax ~ 10~ 3 a) and the same incidence angle (9i — 7r/12). In all three cases, a rapidly increasing ray separation 
is observed, resulting into completely different exit-angles. It will be shown soon that the ray separation increases 
exponentially in time, as the rays bounce back and forth through the dielectric layer. 

In Fig. 3 the exit-angle is shown vs. the (scaled) incidence point position x — x/ a, a being the spatial period of the 
undulated surface, for a fixed incidence angle 6i = n/12, and a few representative values of n/. Intervals of distinctly 
irregular behaviour are observed. The measure of the set of incidence points corresponding to regular behaviour is 
seen to depend on the value of n/, and shrinks to zero as n/ is increased. On the other hand, at fixed n/, the measure 
of the set of incidence points corresponding to regular behaviour is found to be essentially independent of the incidence 
angle, as exemplified in Fig. 4, where nf = 15 and 9i = kir/12, k = 1,3, 5. 
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Closer investigation reveals that smaller and smaller intervals of regular behaviour are hidden within arbitrarily 
small intervals of shown in Fig. 5 for Oi = n/12. 

As anticipated, in the irregular regime, the separation between neighbouring incident rays increases exponentially, 
as the rays bounce back and forth through the dielectric layer. This is shown in Fig. 6, where the logarithm of the 
ray-separation (in units of its initial value) is displayed as a function of time (sec) for U/ = 15, and 6i = kir/12, 
fc = l,3,5. 



III. FULL WAVE ANALYSIS: RESULTS 



A full-wave (Floquet mode) analysis (which is omitted here for brevity, and will be found in a forthcoming paper) 
suggests that in the limit of increasingly large a/Xo, Xo being the field wavelength in vacuum both the internal (z < 0) 
and the external (z > 0) fields eventually comply with Berry's conjecture ||. 

This is illustrated in Fig.s 7 and 8, which refer to the special (scalar) case where the electric field has only a nonzero 
component orthogonal to the xz plane, and all other parameters hold the same values as in Sect. 2. 

The field scattered into the half-space z > is indeed a superposition of (propagating or evanescent) plane waves 
whose wave- vectors span the characteristic (grating-lobe) directions |p"lf 

sin# p = sin#i — n/— -, p<EZ. (3-1) 

with (complex) amplitudes T v E j w here Ei is the (complex) amplitude of the incident plane wave, and 9i is the 
incidence angle. It is seen from ( |3.l| ) that the number of propagating scattered plane waves is oc nja/Xo- Numerical 
simulations further show that upon increasing a/Xo, the distribution of the phase among the pertinent (complex) 
coefficients T p becomes more and more evenly distributed in (0, 2ir). 

On the other hand, as seen from Fig. 7, the agreement between the (averaged normalized) spatial autocor- 
relation function (E[r + x(s/2)]E* [r — x{s/2)\) of the internal (z < 0) field and the Jo(n(z)ko\s\) Bessel function 
limit-form predicted by Berry's conjecture is remarkably good already at a/Xo = 6. 

It is also noted that the complex coefficients T p show a strong sensitivity to the incidence angle 6i, as illustrated 
in Fig. 8 for a/Xo = 9. In Fig. 8 the case of a plain sinusoidal surface is also shown (rif = 1), for which the r p s are 
seen to depend fairly smoothly on the incidence angle. This peculiar behaviour has a special notable consequence, 
also suggested by numerical simulations. In the limit where a/Xo ~ * the spread of the distribution of field values 
at any given point corresponding to a large number of different incidence angles in a sector of aperture A9i. also 
remains 0(A#j), and the mean changes with the chosen field point. On the other hand, in the opposite limit where 
a/Xo ^ 0Oj f° r o,ny nonzero A0j, albeit small, the spread of the distribution blows-up to a finite value, while both 
the spread and the mean of the distribution become almost independent from the chosen field point. This scenario is 
likely to occur in all wave-scattering problems exhibiting chaos in the A — > (ray) limit. 



IV. CONCLUSIONS - APPLICATION POTENTIAL 



Far from being a mere exercise of pure academical interest, the proposed scattering system holds the potential for 
several applications of practical interest. 

The findings presented above suggest that an incident narrow coherent beam of electromagnetic radiation would 
be scattered off the system incoherently and isotropically if a/Ao 3> 1. Thus the proposed system might be used as a 
coating for radar-puzzling applications. 

The multi-hop nature of the ray-paths also suggest that upon addition of suitably small |l2|] dielectric and/or 
backing-surface losses the system might act as a smart radio-frequency absorber, featuring both fairly uniform power 
deposition in the dielectric bulk, and nicely small power reflection-coefficient for a wide interval of incidence angles. 
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Captions to the Figures 



Fig 1 - Sketch of problem's geometry. A perfectly conducting undulated surface with peak-to-peak height A and 
spatial period a is topped by a dielectric layer with maximum thickness h. The refractive index of the layer matches 
that of the vacuum at z = and grows exponentially with depth up to a value rif at z = —h. 



Fig. 2 - Multi-hop ray-paths originating from close-by incidence points display rapidly increasing separations, and 
eventually emerge with widely different exit angles. 



Fig. 3 - Exit angle (deg.) vs. scaled incidence point position x = x/a for fixed incidence angle (0, = n/12) and 
different values of n/: n/ = 1.4 (top), n/ = 30 (mid), n/ = 90 (bottom). 



Fig. 4 - Exit angle (deg.) vs. scaled incidence point position x — x/a for fixed rif and different incidence angles. 
9i = tt/12 (top), 6i = 3tt/12 (mid), 0* = 5tt/12 (bottom). 



Fig. 5 - Exit angle (deg.) vs. scaled incidence point position x = x/a for rif = 15 and 0j = 7r/12. Intervals 
corresponding to regular behaviour are observed at smaller and smaller scales (top to bottom). 



Fig. 6 - Natural logarithm of separation between close- by incident rays as a function of time for rif = 15 and different 
incidence angles: 0, = 7r/12 (left), 0j = 37r/12 (mid), 9i = hir/12 (right). 



Fig. 7 - Full wave analysis with a/Ao = 6. Averaged normalized spatial autocorrelation function of the internal 
(z < 0) field. 



Fig. 8 - Full wave analysis with a/\ = 9. Top: modulus (left) and phase (right) of T p for two close-by incidence 
angles. Bottom: the same for a plain (n/ = 1) undulated surface. 




Fig. 2 - Multi-hop ray-paths originating from close-by incident points. 
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